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The performance of Niobium-based Superconducting Radio Frequency (SRF) particle accelerator cavities 
can be sensitive to localized defects that give rise to quenches at high accelerating gradients. In order to identify 
these material defects on bulk Nb surfaces at their operating frequency and temperature, a wide bandwidth 
microwave microscope with localized and strong RF magnetic fields is developed by integrating a magnetic 
write head into the near-held microwave microscope to enable mapping of the local electrodynamic response in 
the multi-GHz frequency regime at cryogenic temperatures. This magnetic writer demonstrates a localized and 
strong RF magnetic held on bulk Nb surface with B^urface > 10^ mT and sub-micron resolution. By measuring 
the nonlinear response of the superconductor, nonlinearity coming from the nano-scale weak link Josephson 
junctions due to the contaminated surface in the cavity fabrication process is demonstrated. 

PACS numbers: 74.70.-b, 07.79.-v, 74.25.N-, 74.25.nn, 74.50.-l-r 


I. INTRODUCTION 

There is continuing interest to build high energy acceler¬ 
ators to bring electrons and positrons up to near light veloc¬ 
ity using bulk Niobium (Nb) superconducting radio frequency 
(rf) cavities ffl. These electromagnetic cavities resonate at 
microwave frequencies and can efficiently transfer microwave 
jOTwer into the kinetic energy of the charged particle beam 
13113- Each cavity operates with regions of high electric and 
magnetic rf fields, with the highest electric field on the ac¬ 
celerating axis and maximum magnetic field on the equator 
surface of the cavity JH]. However to put these lab-fabricated 
cavities into commercial mass-production, many issues arise. 
It is a challenge to manufacture such high performance cavi¬ 
ties with a high degree of consistency. Much debate regarding 
this inconsistent performance concluded that certain types of 
defects on the Nb cavity surface can behave as a source of 
quenching in tangential high RF magnetic fields 01llllllSl- A 
quench is when the Nb superconductor returns to the normal 
state in the presence of strong fields, thus limiting their util¬ 
ity. Candidates for these defects include pits, oxides, hydrides, 
impurity inclusions, grain boundaries etc. M. These de¬ 
fects arise from cavity fabrication steps such as forming, ma¬ 
chining, electron beam welding, electro polishing (EP) and 
buffered chemical polishing (BCP). These defects are either 
non-superconducting or have lower Tc (for example, NbOj: 
with X = 0.02 ~ 0.04 has Tc ranging from 5.1 K to 6 K i^), 
making them sources of dissipation and interrupting the su¬ 
perconducting current flow. Unfortunately, it is difficult to to¬ 
tally remove all of these defects even after sophisticated phys¬ 
ical and chemical treatments. Because different defects have 
their individual quench limit, not all of the defects behave as 
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sources of quenching in certain operation conditions (for ex¬ 
ample, a 9 cell TESLA cavity usually operates at a temper¬ 
ature of 2 K, an initial accelerating gradient of 25 MV/m, 
1.3 GHz RF frequency under 16.5 MW microwave power 
j^). Unfortunately, it is not known which defects have lower 
quench limits and therefore it is necessary to quantitatively 
understand each type of defect. 

Due to the defect issues with Nb, ideally one would like 
a microscopic technique that identifies defects based on their 
poor microwave performance at low temperatures in the su¬ 
perconductive state. Quenches often occur in regions of the 
cavity with strong tangential magnetic field, hence a local 
probe with strong rf magnetic field concentration is desired. 
Generally speaking, creating a concentrated microwave mag¬ 
netic field in a small area of the sample is a challenge and has 
created interest among several research groups using different 
methods. Fig. [T] shows a brief summary of developments of 
localized high RF magnetic field measurements. Our objec¬ 
tive is to move into the upper-left quadrant of this diagram, 
in other words to develop a microscopic probe that simultane¬ 
ously produces a very strong rf magnetic field over a very lim¬ 
ited area of the sample. The conventional method is to make a 
loop-shaped coil as small as possible to generate a strong mag¬ 
netic field on the sample if lffliTll]ifl^ioifT^ . Some groups 
also bundle many loop-shaped coils to generate a strong lo¬ 
calized field 1131, but unfortunately this method will sacri¬ 
fice spatial resolution. Recently several groups have used 
the method of either photo-lithography or E-beam lithogra¬ 
phy to fabricate a micro-coil and successfully generated about 
10 to 50 Tesla of pulsed magnetic field on 30 ns time scales 
IS 10. One group also demonstrated that this nanosecond 
pulsed field can switch the magnetization of nano magnetic 
particles ifTsIl . but this method will generate heat and may not 
be useful to probe the properties of superconductors due to the 
complicated sample heating during the measurement. Another 
study on utilizing spin-torque-driven microwaves to generate 
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FIG. 1. High-Field, High-Resolution AC/Transient Magnetic Microscopy Generating Highly-Confined Magnetic Fields with Frequency > 
100 MHz : A selected summary of methods to generate a highly concentrated RF magnetic field. Each technique lists a typical RF magnetic 
field scale, degree of spatial confinement/resolution, and an associated time/frequency scale demonstrated. All of the candidate technologies 
are numbered and referenced as follows (1) S. C. Lee ^ loop probe, fldl jn]] llTl . (2) G. Lamura =)> pancake coil d (3) T. Takamasu ^ 
micro-coil d (4) M. Bonfim => micro-coil iia (5)N. Kikuchi =► micro-coil (H (6) G. Woltersdorf => microwave assisted recording d. 
(7) T. Tai => magnetic write head i2^f2§i\ . (8) K. Yoshida =)> spin-torque-driven oscillator (^ . (9) D. Fang =)> spin-torque-driven oscillator 
IHIl . (10) X. Zhu =)> spin-torque-driven oscillator 1^ . Note that larger circle (#) and star sign (■^) indicates the operation frequency/speed 
corresponds to higher frequencies. The star sign also indicates the method used in this paper. 


a localized RF magnetic field to flip the magnetization direc¬ 
tion of a thin single domain mag netic element has been widely 
discussed recently ifl^ll^llzilll^llzsl] . However this tech¬ 
nique doesn’t have application yet in investigating the proper¬ 
ties of SRF cavities. 

In the microwave microscope approach described here, a 
magnetic write head from a hard disk drive is used as a 
near field scanning probe to generate strong and localized 
RF field on a superconducting sample iSl llz^ llzsll . 
The device is a high resolution magnetic microscope probe in 
the GHz frequency regime. The superconducting sample re¬ 
sponds by creating screening currents to maintain the Meiss¬ 
ner state in the material. These currents inevitably produce a 
time-dependent variation in the local value of the superfluid 
density, and this in turn will generate a response at harmonics 
of the driving tone. At localized defective regions, which will 
cause pre-mature quench of the cavity, the enhanced harmonic 
response will be more significant. 

The nonlinear third-harmonic response frequency is studied 
in this paper because that one is intrinsically always present 
(to greater or lesser extent) in a superconductor in zero dc 
field at finite temperature. In addition to this intrinsic re¬ 
sponse, extrinsic properties (i.e. associated with defects) gen¬ 
erate their own third harmonic response, almost always at a 
level far stronger than the intrinsic response 1^ . Figure|2](a) 
illustrates a schematic situation in which defects stimulated by 
a localized scanned RF magnetic field probe show enhanced 
nonlinear response compared to the background intrinsic re¬ 


sponse. One can parameterize the nonlinear response with an 
associated “nonlinear critical current density” called 7jVL- The 
nonlinearity is quantified in terms of the effect that the defect 
has on the superfluid density n, with increasing RF current J. 
This is written as = 1 — ij^)^ , showing that regions 

with small 7jVL will have enhanced suppression of the super¬ 
fluid density ifl^fill]ifl^ll^ . The values of Jnl associated 
with the intrinsic Nb background material, and several defects, 
are also shown schematically in Fig. |2](b). The third harmonic 
power is inversely proportional to Jnl as F*?/ ~ 
film superconductors fl^fllllll^ . Hence defective regions 
with smaller Jnl will produce larger harmonic power (F*?/), 
revealing their presence in the scanned probe microscope. 

This paper will first address experimental measurements 
(section II. and section III.) of the nonlinear third harmonic 
response from bulk Nb materials in localized regions down to 
sub-micron length scales. In section IV, a detailed discussion 
of the quantitative modeling of the experimental results will 
be given and then section V contains the conclusion. This pa¬ 
per provides a more detailed examination of the preliminary 
data presented in Ref. Q. 


II. EXPERIMENT 

A bulk Nb sample with of 9.2 K and well-defined loca¬ 
tions of defects is an ideal candidate for measurement of non- 
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FIG. 2. (a) Schematic Ji^i variation on the superconducting surface 
from intrinsic Nb to different discrete defects. The generation 
from defects is local. The defects with greatest disruption of the 
superfluid density have the lowest and strongest P-^f. Hence 
measuring P 3 f at different positions can be used to find defects on 
superconductors, (b) A schematic hierarchy of a few representative 
defects and their associated values. 



FIG. 3. (a) A schematic picture of the bulk Nb tensile specimen 
before deformation. All dimensions of the specimen are labeled, (b) 
The finished specimen after tensile test and welding treatment in the 
center of the sample. 


linearity using a microwave microscope. Fig. |3] (a) shows a 
schematic diagram of the original tensile specimen made from 
bulk Nb with a high residual-resistivity ratio (RRR). The sur¬ 
face of three tensile specimens were mechanically polished 
and then electropolished prior to deformation. To simulate an 
elliptical SRF cavity forming a welding process using large 
grain niobium, single crystal samples were strained to 40% 
elongation, cut in half, and different halves were welded to¬ 
gether. Fig. [3 (b) shows the finished specimen. The sam¬ 
ple has been carefully imaged by electron backscattered pat¬ 
tern (EBSP) orientation imaging microscopy (OIM) ll^ and 
a variety of microstructural conditions were identified, includ¬ 
ing as-deformed regions that were kept cold near the end, re¬ 
crystallized regions with several grain boundaries in the cen¬ 
ter, and recovered regions in between, many niobium oxides. 



FIG. 4. The setup inside the cryostat probe station for measuring the 
bulk Nb nonlinearity. The slashed region on the top of the bulk Nb 
sample indicates where the P-^f measurements are performed. The 
circled zoom-in picture schematically highlights the magnetic write 
head probe on the surface of the bulk Nb. Note that thermal anchor¬ 
ing is critical for cooling the bulk Nb sample below Tc- 


grain boundaries and dislocations on the surface. The OIM 
analysis identified several specific defects that appear in the 
processing of bulk Nb cavities. 

In order to make the sample surface superconducting in our 
microscope cryostat system, thermally anchoring the sample 
and good radiation shielding from the outside environment is 
crucial. Figure|4]shows the setup inside the cryostat. The bulk 
Nb sample is clamped in place by two halves of an oxygen- 
free high conductivity (OFHC) copper bar. Indium is filled 
within the gap between the Nb sample and OFHC copper 
bar to reduce the thermal resistance between bulk Nb and the 
OFHC copper bar. The two halves of the bar are bolted to¬ 
gether laterally by two 2 — 56 hex-head screws from the side 
and then the whole OFHC copper bar assembly is bolted di¬ 
rectly onto the surface of the sample stage in the cryostat with 
vacuum grease in between to reduce the thermal resistance. 
Note that the sample stage is also made of OFHC copper and 
is also firmly bolted on the top of the cold plate. A thermome¬ 
ter is placed on the bulk Nb top surface to monitor the temper¬ 
ature of the relevant surface of the sample. An OFHC copper 
cryogenic radiation shield is attached to the cold head of the 
cryostat. The lid of the radiation shield has a 3 inch diame¬ 
ter infrared blocking window as a view port. This setup can 
cool down the bulk Nb sample top surface to 5.0 K without 
pumping the exhaust from the cryostat. By pumping the ex¬ 
haust, the surface temperature of the bulk Nb can go below 
4.0 K. However 5.0 K is cold enough to analyze the surface 
microwave properties of this bulk Nb sample. Note that if 
these extensive precautions are not taken, it is virtually impos¬ 
sible to cool the surface of bulk Nb down below its transition 
temperature. More photos inside the vacuum chamber can be 
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FIG. 5. The P-^f dependence on Pf at a selected temperature, 8.7 K, 
for local microwave excitation of bulk Nb with the magnetic write 
head probe. Note the Tc of the bulk Nb sample is 9.2 K. The P^f at 
T =20K results from the probe nonlinearity itself. 


found in the reference | 

The principles of measuring the complex third harmonic 
voltage or related scalar power that is gen¬ 

erated by the sample due to local excitation at a fundamental 
frequency from the magnetic write head has been widely dis¬ 
cussed in previous publications ll^ ll2^ . The mag¬ 
netic write head is made by Seagate (part GTS) and was orig¬ 
inally designed for longitudinal magnetic recording. Hence it 
creates a strong tangential RF magnetic field on the surface 
of the sample. The superconducting sample in proximity to 
the probe responds by creating screening currents to maintain 
the Meissner state in the material. These currents inevitably 
produce a time-dependent variation in the local value of the 
superfluid density, and this in turn will generate a response at 
harmonics of the driving tone. A high pass filter is used out¬ 
side the cryostat to suppress the reflected Pf.if and pick up the 
P^f signal. 


III. THIRD HARMONIC MEASUREMENT RESULTS ON 
BULK NB 


Many localized third harmonic measurements are per¬ 
formed in the slashed region indicated in Fig. 0] Based on 
OIM images, this region has a moderately high density of dis¬ 
organized clusters or knots of dislocations due to the presence 
of some plastic deformation from the tensile test. The probe 
height is determined by optical microscopy from the side of 
the cryostat combined with Sii{f) measurement of the mag¬ 
netic write head to decide the best excitation frequency. The 
probe height is approximately within I jxm from the bulk Nb 
surface. The details of the probe height control are discussed 
in Ref. ll^ll^ . Different locations show slightly different 
microwave nonlinear properties and will be illustrated quali¬ 
tatively in this section. 

Figure |5] shows one of the representative measurements for 


the power dependence of Pt, f with respect to fundamental in¬ 
put power Pf at two fixed temperatures (above Tc and below 
Tc) under 5.025 GHz excitation at this area. The representative 
curve measured atT — 8.7 K (below Tc) shows a distinct P^f 
onset from the noise floor of the spectrum analyzer followed 
by a continuous increase of nonlinearity until a turnover oc¬ 
curs at high excitation power. Then the Pj,f{Pf) data at high 
excitation power eventually approaches and oscillates around 
the probe background nonlinearity, pprobe^Pf)^ measured at 
T = 2QK. Firstly, the onset of the nonlinear response can 
be defined as a temperature dependent lower critical power. 
Pci- This onset power, P^i, will shift toward larger excitation 
power at lower temperature on the Nb superconductor (Fig. 
© lH^. In addition, after a turnover, all P^f{Pf) curves tend 
to decrease with increasing power, suggesting that Nb super¬ 
conductivity is suppressed. The linear extrapolation of P 3 / to 
the noise floor of the spectrum analyzer can be defined as an 
upper critical power, Pc 2 , which suggests that superconduc¬ 
tivity will eventually be destroyed due to the high RF mag¬ 
netic held. It is seen that higher Pc 2 is required to destroy 
the superconductivity at lower temperature 1^ . The tem¬ 
perature dependent Pd and Pc 2 can be associated with tem¬ 
perature dependent surface held Bd and 5^2 (or perhaps Be). 
The analysis of the temperature dependent Pc 2 has been used 
to estimate the surface held excited by the magnetic write 
head probe. This analysis has been carried out by a ht of 
the experimental temperature dependent critical held Bc{T) 
from each experimental Pc 2 with the approximate equation 

Bc{T)^Bc{QK)(j - {T/Tc)^y valid near Tc IH. The ht in¬ 
dicates localized magnetic held from the magnetic write head 
probe is on the order of the thermodynamic critical held of 
Nb with ^surface ^ lO^mT i^ . Different positions will show 
different Pd and Pc 2 values due to the inhomogeneity of the 
bulk Nb surface properties. 

If the inhomogeneity of the Nb surface contains Nb hydride 
or Nb oxide, which has lower critical temperature there may 
form a weak link Josephson junction. In this case, the curve 
of P^f versus Pf will behave different from the case of the 
homogeneous bulk. Fig. | 6 ] shows the P ^/ versus Pf at 4.7 
K comparing to that at T=9 K and T=6.7 K. One can see at 
T=4.7 K, the curve of Pj ,/ versus Pf shows an additional peak 
while the excitation power is between -2.5 dBm and 7 dBm. 
This implies some additional source of nonlinearity is excited 
at this temperature. Detailed interpretation will be given in the 
section IV. 

The temperature dependence of 7 * 3 /(T) at the same probe 
position is shown in Fig. |7] at two values of Pf. The onset 
of Nb nonlinearity begins immediately at the bulk Nb Tc (9.2 
K) and then quickly saturates before T =8 K for Pf=5.5 dBm, 
5.025 GHz excitation. The nonlinearity shows a dip around 
5.8 K. Below 5.8 K, the nonlinearity increases again. This 
same qualitative behavior also happens at a lower excitation 
power of Pf=Q dBm. From this temperature dependent nonlin¬ 
ear behavior, the measured response at high power is consis¬ 
tent with the nonlinearity from the Josephson effect lOj^ . 
This situation will happen on the bulk Nb surface due to the 
presence of oxides and hydrides forming weak link Joseph- 
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FIG. 6 . The P 3 y dependence on Pf at selected temperatures on the 
bulk Nb sample with the magnetic write head probe. The y depen¬ 
dence on Pf at 4.7 K shows a different nonlinear behavior compared 
to the other two temperatures. The nonlinearity of the probe is mea¬ 
sured at T=10.7 K in this data set. 



FIG. 8 . Schematic illustration of an induced screening current loop 
created by the magnetic field pointing out of the paper. Multiple de¬ 
fects such as subgrains are connected by weak link Josephson junc¬ 
tions between the defects. Note 5 is the area of the screening current 
loop. The arrow directions illustrate the direction of diamagnetic 
current. This figure indicates 10 Josephson junctions. 



FIG. 7. The temperature dependent Pif{T) (measured as temperature 
decreases) of the bulk Nb sample at two specific input powers and 
5.025 GHz excitation. Note the Tc of the bulk Nb sample is 9.2 K 
and the onset of nonlinearity at low excitation power also begins at 
exactly 9.2 K. 


son junctions for example. In the next section, the 

nonlinear mechanisms active in the bulk Nb measurement are 
modeled by the weak-link Josephson nonlinear effect. 


IV. MODELING NONLINEARITY 

Possible extrinsic nonlinear mechanisms for the bulk 
Nb measurements includes moviris vortices iIjtI] and the 
weak-link Josephson effect 11^ 11^ . The data presented 
in section III can be reasonably well explained by the 
nonlinear electrodynamics of the bulk Nb superconduc¬ 
tor suffering phase slips ifssi] across a Josephson junction 
lH^]. Josephson junctions will form on the inner sur¬ 
face of the bulk Nb cavities due to the inevitable oxida¬ 


tion in air ll^ and exposure to water vapor. This sur¬ 
face oxidation deteriorates the superconductivity because of 
the formation of superconductor-insulator-superconductor (S- 
I-S), superconductor-normal metal-superconductor (S-N-S) 
or superconductor-constriction-superconductor (S-C-S) struc¬ 
tures i^ . In addition, Josephson junctions may also develop 
on the Nb surface due to the chemical surface treatment during 
fabrication. 

To model the weak-link Josephson nonlinearity, we con¬ 
sider a perpendicular magnetic field component that induces 
a single screening current loop passing through N identical 
defects, which are arrangements of dislocations, dislocation 
tangles or impurities as shown schematically in Fig. These 
sub-grains are weakly connected together in a closed loop of 
area S. This loop can be driven by a combined DC field (B^c) 
and RF field (Bry) on the superconducting surface giving rise 
to a total flux of 


= <Pdc + ^rfsin{m) 

SBdc ^ ^ SBrf ( 1 ) 

with 9dc =- and ^rf = - - 

<I>o 

where <I>(f) is the total normalized magnetic flux through the 
loop, in units of the flux quantum <I>o^ ^dc is the normalized 
DC flux. <I>r 7 and B^f are the normalized flux and amplitude 
of the normal surface RF held, respectively. In this experi¬ 
ment, ^dc is assumed to be zero because there is no applied 
DC held and we assume that there is no stray DC held from 
the write head or the environment. A loop of N defects un¬ 
der the applied RF magnetic flux (Or/) will have an internal 
energy (E) which can be defined by the Josephson coupling 
Hamiltonian in. Clearly, if Ory is higher than the cusp of 
Josephson coupling energy (either higher than the first energy 
cusp, the second or even higher energy cusp), a phase slip will 
occur in the junction . Hence, for the amplitude of RF flux 

shown in Fig. |9] a phase slip will occur at 4 times in the first 
half of the RF cycle, at times t\, t 2 , h and 14 , where fi, 14 , are 
the times when the applied magnetic flux is equivalent to the 
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FIG. 9. (Upper frame) Normalized magnetic flux without the pres¬ 
ence of DC field (<I>(f) = ^rfsin{coi)). This figure assumes the cusps 
of the energy appear when <I> = 0.25 and 0.75. (Bottom frame) The 
generated voltage spikes in a half of RF cycle (t) due to the presence 
of phase slips of Josephson junctions in the loop of N subgrains. 


Z 



FIG. 10. The screening current (J) induced in a superconductor due 
to a point magnetic dipole Tit with a distance “a” above the supercon¬ 
ducting surface plane (Z = a), p and 6 are the radial distance and 
azimuth angle respectively. Surface current calculation from Ref. 

Si. 


first cusp of the energy of the N grain junction ring and f 2 and 
f 3 are the times when the applied magnetic flux is equivalent 
to the second cusp of the energy of the N grains. 

The response voltage due to the phase slip of the Josephson 
junction is given by 


y(f) 


h d<j) 
2 e dt 


( 2 ) 


where <j) is the gauge-invariant phase difference across a 
Josephson junction. This voltage is a series of 4 spikes in 
each half-RF period, corresponding to each of the 4 phase slip 
events. We project out the third harmonic voltage by making 
a sine and cosine expansion in harmonics of the fundamental 
drive frequency. 


2 f'^ 

V$fa = — / V(t)sm(3cot)dt 

T Jo 

2 

Vjfi, = — / V(t)cos(3cot)dt 

T Jo 


(3) 


where T is the rf period. This yields the following expression 
for the 3'^^ harmonic voltage for the case of n times phase-slip 
events in each of period: 


Vifa 

V3fb 


Pi 27Z(0 1 I . /_ , I 

Lk*«(3®h)l 

e ly L . 

Ti '2.71(0 rT—t , 


(4) 


where N is the number of equivalent Josephson junctions in 
the loop. 

In the experiment, the current loop couples flux back to the 
magnetic write head probe. Hence what we measure is an 
induced voltage in the magnetic write head probe given by 
Ve = d^e/dt where <I>g is the time dependent induced flux at 
the write head probe. To estimate the voltage induced on the 


magnetic probe, one can consider a magnetic dipole moment 
over the superconductor at a distance “a” from the supercon¬ 
ducting surface lf3^ iifll] llddl] as shown in Fig. [TO] The 

magnetic dipole can be described by a Dirac delta function 
M5(.r)5(y)5(z — a), where M is the moment of the dipole. 
The superconductor maintains itself in the Meissner state and 
generates a screening current on the superconducting surface. 
The radial and azimuth angle components of the generated 
screening current due to this dipole can be expressed as 14^ 


JeiP.e) 


Jp{p,0) = 


MsinO yjp^ + a^ — a 


McosO + a^ — a a ^ 

47rA2 ^ ~ (p2 + a2)3/2J 


(5) 


where p and 9 indicate the radial distance and the azimuth 
angle with respect to an origin on the surface just below the 
dipole (see Fig. [loll. Hence, the maximum current, Jmax, in¬ 
duced by a magnetic dipole moment M at a height “a” above 
the superconductor with penetration depth A will be at the lo¬ 
cation of p = 0 and 9 — —njl. 


Jmax — Jp{p — 0,9 — n/2) — 

Jg{p=0,9 = -7t/2)=0 


(6) 


The magnetic dipole moment M also can be expressed in 
terms of the magnetic flux in the gap of the probe as M = 
<Pelgap/P 07 where Igap is the length of the magnetic gap on the 
bottom of the magnetic yoke 1231. Hence the induced voltage 
on the probe can be written as 

^ _ d^e d poM dJ STzpoX^a^ 

dt dt Igap dt ^gap 

Take an estimate of the 3'^'^ harmonic current density to be 
J 3 =l 3 l(X a), where I 3 is the third order harmonic current, and 
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in combination with Eq|4] the induced third harmonic volt¬ 
age for sine and cosine expansion and the corresponding third 
harmonic power generated by this phase slip mechanism in 
the weak-link Josephson junction is given in the frequency 
domain by 


nfa = 


96iiQ(oX{T)ah{T) 


Nl 


gap 




96^coX{T)ah{T) 


Nl 


gap 


^ I sin{3(0ti) I 

i 

n 

^ I cos{30)ti) I 


(8) 




H.) +('' 3 '/.)' 

2Zo 



where Zo=50 the characteristic impedance of the mi¬ 
crowave generator. A complete calculation of the tempera¬ 
ture and power dependent has been carried out in refer¬ 
ence im. The outcome shows that I 3 = Ic{T)J^{f3)sin{3(£)t) 
with j3 = coming from the third order expansion of 

the Fourier series for I{t) = Icsin{^sin{3(£)t)) in terms of the 
Bessel function y„(j3) Here 4 is the junction critical 

current, which we take to be the Ambegaokar and Baratoff 
estimation ll^ . 


Ic[T) = 


nA{T) 

2eRn 


tanh[ 


A(r) 

2kBT' 


(9) 


where A(7’) is the temperature dependent superconducting 
gap parameter, and is the normal resistance of the junc¬ 
tion. The temperature dependence of A(r) can be calculated 
according to BCS theory with the zero temperature London 
penetration depth of Nb Aq = 40 nm H. 

To simplify the calculation for fitting experimental P 3 f{Pf) 
and P-if{T) data, the number of Josephson junctions N as 
shown in Fig. 0 in the current loop is assumed to be ten in 
Eq. This assumption is based on Halbritter’s estimation 
that weak links on the surface of bulk Nb are about 10 ^ 100 
nm apart . In addition, both square bracketed terms in 

Eq. [ 8 ] are assumed to be on the order of one. The distance 
between the probe and sample, a, is assumed to be 100 nm, 
which is the scale of the magnetic gap and flying height dur¬ 
ing the reading and writing process in high speed magnetic 
recording Si m. Two parameters, R„ and j}, are left for 
fitting. 

From the power dependence measurement shown in Fig. |5] 
and Fig. | 6 ] one can see that the probe nonlinearity, 
participates in the measurement at high excitation power. In 
order to extract the response from only the superconducting 
sample at high power, an effective F 3 / is defined by taking the 
absolute value of the difference between the measured F 3 / in 
the superconducting state and the measured f’ 3 f in the normal 
state, without considering the relative phases of the contribut¬ 
ing nonlinearity. This crude subtraction process provides a 
qualitative measure of the high-power nonlinear response of 
the sample. Therefore, the effective F 3 f can be regarded as the 
nonlinearity arising from just the superconducting sample, to 
a first approximation. 


FIG. 11. Effective P^f versus Pf for the bulk Nb sample measured 
by the magnetic write head probe at 9 K under 5.025 GHz excitation 
as shown by the dot points. The solid line is the theoretical calcula¬ 
tion based on the phase slip model of the Josephson junction. The 
parameters are as follows: T=9.0 K, Frequency=5.025 GHz, R„ = 90 
5 = 9.11 * lO^*'* frP', a = 100 nm, N=10, Igap = 100 nm. Only /?„ 
and S were varied to perform the fit. 


The points in Fig. [TT]show the effective P^ f at the represen¬ 
tative temperature of 9 K from the experimental data in Fig. 
| 6 ] The solid line in Fig. [TT] indicates the fit to the phase slip 
model for the temperature of T=9.0 K. The relation between 
the incident fundamental power, Pf, and the actual surface 
magnetic field experienced on the bulk Nb surface is taken 
as Pf — kBrp with k — 25.6W directly for fitting the 

effective 4 ?/, where Brf is the amplitude of the surface mag¬ 
netic field on the bulk Nb superconductor and k is an coupling 
coefficient between the probe and sample to connect power 
and magnetic field. The value of k can be judged by the field 
scale generated by finite element simulator and experimental 
results on a known conventional superconductor for calibra¬ 
tion d^. Again, the fits assume constant values of the probe 
height a = 100 nm, the number of Josephson junctions A = 10, 
and the magnetic gap of the yoke Igap = 100 nm. Then only 
the product of 5 = 1 = B^f (which comes from the definition of 

j3 = and /?„ will affect the calculated P^f. 

Finally, two parameters, R„ and S, are varied to fit the model 
[Eq.[ 8 l to the data. Varying/?„ can adjust the amplitude of f’ 3 / 
and varying S can laterally shift the theoretical curve to locate 
the first P$ f {Pf) dip position in the experimental data. From 
Fig. [n] the dip position of the model at = 11.3 dBm and 
Pf —15.0 dBm are matched to dip position of the effective P 3 / 
measurement taken at 9 K under 5.025 GHz excitation by tak¬ 
ing S = {320 nm)^. However, the model shows deep dips but 
the experimental data does not. It is clear in Fig| 6 ]that the data 
of P^f at T < Tc intersects with the Py°^^{T = 10.7K) (back¬ 
ground level) at high excitation power, implying deep minima 
in the harmonic response of the superconductor. The crude 
background subtraction method from limited data points ob¬ 
scures these minima. The amplitude of P 3 f at the first peak 
is matched to the model by taking R„ = 90 H. This junction 
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FIG. 12. Temperature dependent fitting parameters Rn (normal resis¬ 
tance of the junction) and 5 (area of the induced screening current 
loop). Each value is determined hy the fitting method of the phase 
slip model from the P 3 y versus Pf data in Fig. Note that all data 
points for both R„ and S are connected hy an interpolation spline 
function. 


resistance is consistent with Halbritter’s estimation for a bulk 
Nb weak link junction, which is R^i 10^^^ Q.np- Hsl] JUt] . 
At higher power, the amplitude of the effective y does not 
match the phase slip theoretical curve after the first dip. Sev¬ 
eral issues may be at play ; first, when the applied field is 
higher than the vortex nucleation field, generation, motion and 
annihilation of vortices will contribute a significant nonlinear 
response H^l- In addition, vortices pinned by defects will os¬ 
cillate with the competition of pinnin g an d image attraction 
forces, generating harmonic response 11^ . This vortex pin¬ 
ning phenomenon becomes significant at strong fields which 
push the vortex deeper past pinning sites. Another reason is 
that the probe nonlinearity is over-simplified in the calcula¬ 
tion: the phase of the is not included when subtracting 

the measured y in the normal state from the measured P ^y in 
the superconducting state. 

The loop area S for the fit at T=9 K (Fig. [TT]) is 5 = 
9.12* 10^*^ nP', equivalent to a circular loop with radius 
r = 170 nm. This value of the fit loop radius is on the scale 
of the magnetic gap ( ^ 100 nm) in the write head probe and 
implies the resolution of the magnetic write head probe in the 
near field microwave microscope. This length scale is also 
consistent with those put forward by Halbritter for the dis¬ 
tance between weak links on the surface of air-exposed bulk 

Nb iHililiilii. 

A summary of the temperature dependent /?„ and S from 
our model fits for the effective P 3 y at different temperatures 
is shown in Fig. [12] We see that these two parameters are 
temperature dependent, (note that the number of junctions in 
the loop, N, could also vary, but here we choose to keep this 
value fixed). One reason for temperature dependent and S 
is because the screening current loop will look for the easiest 
trajectory with highest critical current of the weak link Joseph- 
son junctions to pass through. One can see the /?„ gradually 


FIG. 13. Temperature dependence of the third harmonic response for 
the hulk Nh sample at the excitation power 5.5 dBm and excitation 
frequency 5.025 GHz measured hy the magnetic write head probe 
(dot points). The noise floor of the spectrum analyzer at this excita¬ 
tion frequency is about -161 dBm. A calculated result (solid curve 
in black) based on the phase slip model is fit to the data. The fitting 
parameter is S^y = 11.75 niT, with fixed values for the remaining 
parameters: lgap=W0 nm, N = 10, a = 100 nm. 


goes up with decreasing temperature from 9 K to 5.7 K. This 
is consistent with the decrease of quasi-particle density upon 
going to lower temperature. For 5(7’), a feature of the jump 
from low-values at T=6.3 K to the larger values at T=5.7 K in¬ 
dicates the loops get bigger at low temperature. The tempera¬ 
ture of this jump is coincident with the transition temperature 
of NbOj: material, or proximity-coupled Nb, in the vicinity of 
weak links in bulk Nb with reduced transition temperatures of 
5.1 K to 6.0 K as identified by Halbritter (|^. Below 5.7 K, the 
Rn{T) decreases with decreasing temperature. This decrease 
also corresponds to a decrease of the area S with decreasing 
temperature below 5.7 K. Although we don’t know how to 
interpret this decrease of /?„ and S at low temperature, sev¬ 
eral possible issues may be in play. The localized heating of 
the sample by the probe may be another possible reason for 
the non-monotonic dependence in Fig. [T2| These possibili¬ 
ties may affect the screening current area S at low temperature 
and result in a corresponding change of /?„ at low tempera¬ 
ture. These results are consistent with the existence of a dense 
network of weak links on the surface of air-exposed Nb. It is 
clear that our microscope is probing this complex network on 
the scale of the defect size and that the network evolves signif¬ 
icantly under RF-field excitation as a function of temperature. 

Finally, the experimental temperature dependent PjfiT) 
with 5.5 dBm excitation in Fig. [Tjis also fitted by the phase 
slip model. An interpolation function for the temperature de¬ 
pendence R„{T) and 5(7’) plotted as the solid line of Fig. [T2| 
and used for P-if{T) fitting. The black solid line in Fig. fOlis 
the fit result from the phase slip model with Brf — 11.75 mT. 
Other parameters (a = 100 nm, N=10, Igap = 100 nm) remain 
the same as the fitting of P^f on Pf. The experimental data 
points are shown in red circles taken at 5.5 dBm excitation 
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power and 5.025 GHz excitation frequency from Fig. |7] Due 
to the finite number of data points in the interpolation function 
for the temperature dependence Rn{T) and S{T), the htting 
curve shows some discrepancy. However, most of the P-ifiJ) 
data points still follow the trend of the phase slip model. The 
ht of the 0 dBm data in Fig. |7]is qualitatively correct as well, 
which makes a compelling case that the nonlinear mechanism 
comes from the weak link Josephson junction while the exci¬ 
tation power is less than the hrst deep oscillation in the effec¬ 
tive P^f{Pf) curve. 

V. CONCLUSIONS 

A clear reproducible nonlinear response signal from the 
surface of superconducting bulk Nb is obtained by the mag¬ 
netic write head probe. The success of nonlinear excitation 
on bulk Nb also implies the magnetic held from the magnetic 
write head probe is on the order of the thermodynamic crit¬ 
ical held of the measured bulk Nb sample. From hts to the 
data, nonlinearity from the weak-link Josephson effect can ex¬ 


plain the measured nonlinear response on the bulk Nb for both 
P-if{T) and Pjf{Pf) measurements. The measured nonlinear¬ 
ity of the bulk Nb samples is perhaps due to the contaminated 
surface in the fabrication process or subsequent handling in 
air. The nonlinear near-held magnetic held microwave micro¬ 
scope has great potential to identify the electrodynamic prop¬ 
erties of bulk Nb materials and to image the electrodynamic 
defects on superconducting Nb in the GHz frequency region. 
Our work in measuring local high frequency microwave prop¬ 
erties of Nb superconductors will beneht the SRF community 
to understand the key issues that lead to degradation in the 
performance of Nb cavities. 


VI. ACKNOWLEDGEMENT 

This work is supported by the US Department of Energy 
/ High Energy Physics through grant # de-sc0012036 and 
CNAM. The work at MSU was supported by DOE-OHEP, 
Contract No. DES0004222. 


[1] Abdelhak Djouadi, Joseph Lykken, Klaus Mdnig, Yasuhiro 
Okada, Mark Oreglia, Satoru Yamashita, et al, “International 
Linear Collider Reference Design Report Volume 2: PHYSICS 
AT THE ILC,” arXiv: 0709.1893 

[2] H. A. Schwettman, J.P. Tumeaure and W.M. Fairbank, T.I 
Smith, M.S. McAshan, F.B. Wilson, E.E Chamber, “Low tem¬ 
perature aspects of a cryogenic accelerator,” IEEE Transaction 
on Nuclear Science 14, (3), 336, PAC (1967). 

[3] Hasan Padamsee, Jens Knobloch, Tomas Hays “RE supercon¬ 
ductivities on accelerators (Wiley Series in Beam Physics and 
Accelerator Technology), ” ISBN: 978-3-527-40842-9, (1998). 

[4] G. Ciovati, G. Myneni, F. Stevie, P. Maheshwari, D. Griffis, 
“High held Q slope and the baking effect: Review of recent ex¬ 
perimental results and new data on Nb heat treatments,” Phys. 
Rev. Special Topics-Accerlerators and Beams 13, 022002, 
( 2010 ). 

[5] Z. H. Sung, A. A. Polyanskii, P. J. Lee, A. Gurevich, D. C. 
Larbalestier, “Suppressed Superconductivity on the Surface of 
Superconducting RF Quality Niobium for Particle Accelerating 
Cavities,” AIP Conf. Proc., 1352, 142-150 (2011) 

[6] A. A. Polyanskii, P. J. Lee, A. Gurevich, Zu-Hawn Sung, D. 
C. Larbalestier, “Magneto-Optical Study High-Purity Niobium 
for Superconducting RF Application,” AIP Conf. Proc., 1352, 
186-202 (2011) 

[7] C Antoine, A Aspart, S Regnault, A Chincarini, “Surface Stud¬ 
ies: Method of Analysis and Results,” Proceeding of the 10th 
Workshop on RF Superconductivity, Tsukuba, Japan (2001). 
CZ Antoine, A Aspart, M Berthelot, Y Gasser, JP Poupeau, 
F Valin, “Morphological and Chemical studies of Nb Samples 
after Various Surface Treatment,” Proceeding of the 9th Work¬ 
shop on RF Superconductivity, 295-303, (1999). 

[8] Z. H Sung, P. J. Lee, D. C. Larbalestier, “ Observation of the 
Microstructure of Grain Boundary Oxides in Superconducting 
RF-Quality Niobium With High-Resolution TEM (Transmis¬ 
sion Electron Microscope),” IEEE Trans. Appl. Supercond., 24, 
68-73 (2014). 


[9] J. Halbritter, “On Extrinsic - weak link - effects in the sur¬ 
face impedance of cuprate - and classical - superconductors,” 
Fifth Workshop on RF Superconductivity, DESY, Hamburg, 
Germany (1991). 

[10] S.-C. Lee, S. M. Anlage, “Study of local nonlinear properties 
using a near-held microwave microscope,” IEEE Trans. Appl. 
Supercond. 13, 3594, (2003). 

[11] S.-C. Lee, S. M. Anlage, “Spatially resolved nonlinearity mea¬ 
surements of YBa 2 Cu 307 bi-crystal grain boundaries,” Appl. 
Phys. Lett. 82, 1893, (2003). 

[12] S.-C. Lee, S.-Y. Lee, S. M. Anlage, “Microwave nonlinearities 
of an isolated long YBa 2 Cu 307_5 bicrystal grain boundary,” 
Phys. Rev. B 72, 024527, (2005). 

[13] D. I. Mircea and T. W. Clinton, “Near-held microwave probe for 
local ferromagnetic resonance characterization,” Appl. Phys. 
Lett. 90, 142504, (2007). 

[14] T. W. Clinton, N. Benatmane,J. Hohlfeld, Erol Girt, “Compar¬ 
ison of a near-held ferromagnetic resonance probe with pump- 
probe characterization of CoCrPt media,” J. Appl. Phys. 103, 
(7), 07F546, (2008). 

[15] G. Lamura, M. Aurino, A. Andreone, J.-C. Villegier “First crit¬ 
ical held measurements of superconducting hlms by third har¬ 
monic analysis,” J. Appl. Phys, 106, 053903, (2009). 

[16] T Takamasu, K Sato, Y Imanaka and K Takehana, “Fabrication 
of a micro-coil pulsed magnet system and its application for 
solid state physics,” Journal of Physics: Conference Series 51, 
591, (2006). 

[17] K. Mackay, M. Bonhm, D. Givord, and A. Fontaine, “50 T 
pulsed magnetic helds in microcoils,” J. Appl. Phys. 87, 1996, 
( 2000 ). 

[18] N. Kikuchi, S. Okamoto, O. Kitakami, “Generation of nanosec¬ 
ond magnetic pulse held for switching experiments on a single 
Co/Pt nanodot,” J. Appl. Phys. 105, 07D506, (2009). 

[19] Georg Woltersdorf, Christian H. Back, “Microwave assisted 
switching of single domain NimFeon elements,” Phys. Rev. 
Lett. 99, 227207, (2007). 



10 


[20] K. Yoshida., E. Uda, N. Udagawa, and Y. Kanai, “Investiga¬ 
tion on magnetic fields from field-generating layer in MAMR 
IEEE Trans. Mag. 44, 3408, (2008). 

[21] D. Fang, H. Kurebayashi, J. Wunderlich, K. Vyorny, L. P. 
Zdrbo, R. P. Campion, A. Casiraghi, B. L. Gallagher, T. Jung- 
wirth, A. J. Ferguson 1, “Spin-oebir Driven Ferromagnetic Res¬ 
onance,” Nature Nanotechnology, 22, 413, (2011). 

[22] X. Zhu, J. G. Zhu, “Bias-field-free microwave oscillator driven 
by perpendicularly polarized spin current,” IEEE Trans. Mag. 
42, 2670, (2006). 

[23] C. H. Sim, M. Moneck, T. Liew, J.G. Zhu, “Frequency- 
tunable perpendicular spin torque oscillator” J. Appl. Phys. Ill, 
07C914, (2012). 

[24] Tamin Tai, X. X. Xi, C. G. Zhuang, D. I. Mircea, S. M. Anlage, 
“Nonlinear near-held microwave microscope for RF defect lo¬ 
calization in superconductors,” IEEE Trans. Appl. Supercond. 
21, 2615, (2011). 

[25] Tamin Tai, B. G. Ghamsari, T. Tan, C. G. Zhuang, X. X. Xi, 
Steven M. Anlage, “MgB2 nonlinear properties investigated un¬ 
der localized high rf magnetic held excitation,” Phys. Rev. ST 
Accel. Beams 15, 122002, (2012). 

[26] Tamin Tai, B. G. Ghamsari, Steven M. Anlage, “Nanoscale 
electrodynamic response of Nb superconductors,” IEEE Trans. 
Appl. Supercond. 23, 7100104, (2013). 

[27] Tamin Tai, B. G. Ghamsari, Steven M. Anlage, “Modeling the 
nanoscale linear response of superconducting thin hlms mea¬ 
sured by a scanning probe microwave microscope,” J. Appl. 
Phys. 115, 203908 (2014). 

[28] Tamin Tai, Behnood G. Ghamsari, Thomas R. Bieler, Teng Tan, 
X. X. Xi and Steven M. Anlage, “Near-held microwave mag¬ 
netic nanoscopy of superconducting radio frequency cavity ma¬ 
terials,” Appl. Phys. Lett. 104, 232603 (2014). 

[29] D. E. Oates, Y. D. Agassi, B. H. Moeckly, “Intermodulation 
distortion and nonlinearity in MgB 2 ' experiment and theory,” 
IEEE Trans. Appl. Supercond. 17, (2), 2871-2874, (2007). 

[30] J. C. Booth, J. A. Beall, D. A. Rudman, L. R. Vale, R. H. Ono, 
“Geometry dependence of nonlinear effects in high temperature 
superconducting transmission lines at microwave frequencies,” 
J. Appl. Phys. 86, 1020, (1999). 

[31] Di Kang, Derek C. Baars, Thomas R. Bieler, and Chris C. 
Compton, “Characterization of large grain Nb ingot microstruc¬ 
ture using EBSP mapping and Laue camera methods,” AIP 
Conf Proc. 1352, 90-99, (2011). 

[32] Tamin Tai, “Measuring electromagnetic properties of supercon¬ 
ductors in high and localized RF magnetic held,” Ph.D. dis¬ 
sertation, University of Maryland-College Park, (2013), see 
http://hdl.handle.net/1903/14668 

[33] Carson Jeffries, Q. Harry Lam, Youngtae Kim, L. C. Bourne, 
and A. Zettl, “Symmetry breaking and nonlinear electrodynam¬ 
ics in the ceramic superconductor YBa2Cuj,0-]T Phys. Rev. B. 
37, 9840, (1988). 

[34] Ting-kang Xia, D. Stroud, “Nonlinear electrodynamics and 
nonresonsnt microwave absoption in ceramic superconductors,” 


Phys. Rev. B. 39, (7), 4792, (1989). 

[35] J. Halbitter, “On the oxidation and on the superconductivity of 
niobium,” Appl. Phys. A-Mater. 43, (1), 1-28 (1987). 

[36] G. Ciovati and J. Halbritter, “Analysis of the medium held Q- 
slope in superconducting cavities made of bulk niobium,” Phys- 
ica C: Superconductivity 441, (1-2), 57-61 (2006). 

[37] A. Gurevich, G. Ciovati, “Dynamics of vortex penetration, 
jumpwise instabilities, and nonlinear surface resistance of type- 
II superconductors in strong rf helds,” Phys. Rev. B 77, 104501, 
(2008). 

[38] Michael Tinkham, “Introduction to superconductivity,” 2nd 
Edition, ISBN-10: 0486435032 — ISBN-13:978-0-486-43503- 
9, (2004). 

[39] M. V. Milosevic, F. M. Peelers, “Vortex pinning in a supercon¬ 
ducting him due to in-plane magnetized ferromagnets of dif¬ 
ferent shapes: The London approximation,” Phys. Rev. B. 69, 
104522, (2004). 

[40] A. V. Kapra, V. R. Misko, D. Y. Vodolazov and F. M. Peelers, 
“The guidance of vortex-antivortex pairs by in-plane magnetic 
dipoles in a superconducting hnite-size him,” Supercond. Sci. 
Technol. 24, 024014, (2011). 

[41] M. J. Van Bael, J. Bekaert, K. Temst, L. Van Look, M. V. 
Moshchalkov, Y. Bruynseraede, G. D. Howells, A. N. Grig- 
orenko, S. J. Bending “Local observation of held polarity de¬ 
pendent flux pinning by magnetic dipoles,” Phys. Rev. Lett. 86, 

(I) , 155-158, (2001) 

[42] M. V. Milosevic, S. V. Yampolskii, and F. M. Peelers, “Mag¬ 
netic pinning of vortices in a superconducting him: The an¬ 
tivortex magnetic dipole interaction energy in the London ap¬ 
proximation,” Phys. Rev. B 66, 174519, (2002). 

[43] V. Ambegaokar and Baratoff, “Tunneling between supercon¬ 
ductors,” Phys. Rev. Lett. 10, 486, (1963) 

[44] J. R. Waldram, “Superconductivity of metals and cuprates,” 
Chapter 10. (The non-local form of BCS theory) ISBN-10: 
0852743378 — ISBN-13: 978-0852743379, (1996). 

[45] J. Halbritter, ’’Materials science and surface impedance Z 
(T,f,H) of Nb and YBCO and their quantitative modelling by 
the leakage current of weak links,” Supercond. Sci. Tech. 12, 

(II) , 883-886, (1999). 

[46] Klaas B. Klaassen, Jack C.L. van Peppen, “Nanosecond and 
sub-nanosecond writing experiments,” IEEE Trans. Magn. 35, 
No.2, 625, (1999). 

[47] F. H. Liu, S. Shi, J. Wang, Y. Chen, K. Stoev, L. Leal, R. Saha, 
H.C. Tong, S. Dey, M. Nojaba, “Magnetic recording at a data 
rate of one gigabit per second,” IEEE Trans. Appl. 37, (2), 613- 
618, (2001). 

[48] J. Halbritter, “Granular superconductors and their intrinsic and 
extrinsic surface impedance,” Journal of Superconductivity 8, 
(6), 691-703, (1995). 


